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THERMAL  CALCULATION  OF  A  THERMOGENERATOR  AT  CHANGING 
TEMPERATURES  ALONG  THEHMOCONTACT  S'JR'.-AOES 

G.  A.  Varshavskiy,  I.  A.  F.ezgol' 

( Moscow) 

The  arlicle  finds  expressions  fo”  the  output,  power  of  efficiency 
of  a  thermoe]  ectr  ic  gerjerator  ar.d  tern' crature  distribution  along  the 
heat  earner  under  the  condition  that  ’  cmperature?  of  the  thermocoii- 
tacr,  surfaces  vary  due  to  the  cool  oig  and  heating  of  the  heat  car¬ 
riers.  Simple  approximation  oalouia'.  on  formulas  are  given  for  the 
particular  case.s  examir:ed. 

/ 

In  connection  with  auccGose.s  !'■  semiconductor  technology,  there 
i  .s  widespread  use  ol  ther  moelec 1 1  .1  c  ..onvcrters  for  obtaining  elec¬ 
trical  energy  and  i'or  cooling  purport;?,.  The  therruoteclin i ca  1  calcu¬ 
lation  and  1  invest .1  gat  1  ori  of  the  coo.' ation  of  the  converters  are  usual¬ 
ly  done  by  means  of  elementary  r.hac  etical  analysis  proposed  by  A.  E. 
Ioffe  in  his  time  [I] -and  refined  somewhat  later  by  A.  I.  Burshteyn 
[2]  and  accepted  in  all  the  world  literature  [3-5]. 

The  estate  ished  mcthoci.t  of  calculation  make  it  possible  to  deter¬ 
mine  the  efficiency  of  these  device.s  and  f.i.nd  the  optimal  conditions 
of  tii.rir  operation  at  constant  t.cinperatures  of  hot  and  cold  soloenrig 
This  simplification,  being  scheme  plan,  can  be  justified  by  two  cori- 
.siderations  of  a  practical  nature.  The  first  rcfiecL-s  the  technical 
Side  of  the  question,  since  models  of  thermoelectric  generator’s  (TEG), 
created  at  the  initial  stage  of  development  of  this  branch  of  tech¬ 
nology,  functior.ied  at  practically  constant  terrip‘=' natures  of  hot  and 
cold  jurjctions.  These  are,  for'  example,  the  thermoelectric  gcoiera- 
tor  on  a  ke-usene  lamp  [1-31,  the  thermogenorator  of  appai'atus 
SNA!’-JI1  with  an  isotopic  heat  source  [6],  diffcrenu  solar  T.hcririogcMi- 
erators  [7,  8],  and  .sir  on.  The-  second  circumstance  i  .s  cuiinectcd  with 
the  fact  tiiat  i iicrnioe  1  cc tr  1  c  parameter's  of  tlie  substarjce  greatly 
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depend  on  i.emperanure ,  and  a  detailed  account  of  this  dependence  for 
varial)le  temperature  limits  leads  to  expression  which  cannot  be  in¬ 
vestigated  analytically  to  the  end. 

The  striving  to  increase  the  power  of  the  power  installations 
inevitably  led  to  the  creation  of  a  two-circuit  sc''eme  of  the  thermo- 
generator  19],  where  the  transporting  of  the  heat  is  produced  by  heat 
earners,  the  temperature  of  which  considerably  varies  as  a  result  of 
the  heat  transfer  from  a  hot  liquid  to  a  cold  one  along  the  channel 
of  the  heat  carrier.  With  the  thorough  heat-engineering  calculation 
of  such  apparatuses,  the  dependence  of  thermoelecf^ic  characteristics 
of  the  materials  on  temperature  must  be  considered. 

However,  in  the  preliminary  analysis  of  the  operation  of  the 
uhermogenerator ,  it  is  expedient  to  determine  its  output  parameters 
Without  taking  into  account  the  effect  of  variable  temperatures  of 
the  .solderings  on  thermoelectric  properties  of  the  substances  but 
considei'ing  the  change  in  temperatures  of  heat  carriers  along  the 
t her moGi.li tact  surfaces  . 

This  work  is  devoted  to  the  finding  of  analytical  bonds,  which 
make  it  possible  to  make  tiie  calculation  of  the  ther'mogenerator  at 
changing  (as  a  result  of  heat  transfer)  temperatures  of  the  solder- 
ings.  The  obtained  expressioris  cari  be  useful  iri  the  preliminary  de- 
Lerrnination  of  the  optimal  parameters  of  the  generators  (maximal 
power'  and  efficiency,  mrriimal  weight,  etc.)  and  an  analysis  of  pro¬ 
cesses  of  regulation.  '  •  ■  ■  ''  '  •  '  ■■  ' 

Given  at  the  end  of  the  article  are  certain  recommendations 
according  to  an  analysis  of  the  TKG  operation,  taking  into  account 
the  dependence  ol  ther'nioelectr ic  properties  of  materials  on  tempera- 
t  u  r  ■  e  . 

Calculation  of  characteristics  of  the  thermoelectric  gc’r,erator. 
Let  us  examine  the  THG  operation  in  the  two-cii’cuit  scheme  {.  ig.  1). 
The  hot  heat  carrier  1,  which  circulates  in  the  first  circuit,  feeds 
heat  Ll. rough  the  " thermocontact  packet"  2  to  the-  hot  junctions  [sol- 
derir;gs]  of  the  thermoelements.  I’ar-t  of  thi.s  heat  is  converted  into 
electrical  eticrgy,  and  the  I'Cina  i  n  i  ng  heat,  is  dr-awii  off  fr'om  the  cold 
juiintjuiis  of  ttie  tlicrmeel  ements  through  the  "  I  her-rnocontact  packet" 
to  tt]i_  Cold  h'  itt  cii'i'iCr  b  c  i  rci>  1  t ..  t.g  in  the  se''onn  cjr.juji  .  'li'- 
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Lhst'moeeiierator  ii.'-ludes  a  grear.  number  of  single  thermoel  emeiils , 
which  are  contiecLed  electrically  to  each  other  in  series  or  in  paral¬ 
lel,  depending  on  the  required  output  parameters  of  the  TEG,  The 
ba.sic  dimensions  of  the  thermoelements,  height  h  and  cross  sections 
Sp  and  Sji  of  thi-^  br  ariohes,  are  identical  for  the  whole  thermogenerator . 

In  the  course  of  the  analysis  we  will  use  the  average  valuer-  of 
the  charac ter- 1  St  1 cs  for  the  positive  and  negative  branches  of  the 
thermoelement  (symbols  are  given  at  the  end  of  the  article) 

1  .N,, 

Cl  0(  p  -I  <(,,  )  .  X  --  X,.  t  x„ 

2  M, 

1  '  .S„  \  b'j,  1  „ 

i  V  .Sp  .s.,  *50 

Let  us  take  the  following  simplifying  assump t -  . ns : 

1)  the  electrophy.sical  character sties  n  x,  (i  .  averaged  over 
the  op-.-rating  temperature  range,  are  identical  for  a.l  the  thermo- 

c 1 emen  ts ; 

2)  the  average  current  density  j  is  constant  for  ar,y  cross  sec¬ 
tion  of  the  thermogenerator ; 

3)  the  thermal  i-esistances  and  kj;  of  the  '' t  h'-rnio:  ont  act 
packets"  are  identical  for  ail  the  thermoelements,  and  tne  c-c;rres- 
pondirig  temperature  drops  ai  e  proportional  to  thi  thermtl  fl  jws. 


fr 
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1'  1  g  .  1  . 

A.s  usual  [  1-3],  we  do  not  consider  the  effect  of  the  Thomson 
effect  (since  tlie  specific  thermal  errif  cx  is  aver-aged  over  the  oper¬ 
ating  tomperaturc  range),  and  we  assume  tnat  the  Joule  heat  liberated 
in  the  thennuc  Icmerit  is  transported  equally  to  tlic  liot  and  cold 
Junctions. 
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In  this  case,  for  the  "hot”  and  "cold"  thcrmocontact  surfaces 
of  the  generator  select,ed  along  the  channel  of  heat  carriers  of  ele¬ 
ments  dF  (Fig.  2),  we  obtain  correspondingly  the  following  relations 
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The  solution  to  these  equations  makes  it  possible  to  find  the 
temperature  distribution  along  the  thermocontact  surfaces  of  the  gen 
eratoi',  tiie  fed  and  r'cmoved  quaritities  of  heat,  electrical  power, 
efficiency,  and  so  oti. 
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Fig.  2 . 

H-'--  and  further,  when  tiicn-n  is  a  dual  sign  (  d.  or  ip  )  ,  the 
uppi  '  .sif,  refei^.s  to  the  case  of  cocurrent  flow  aiid  the  lower  sign, 
lo  t'n-  caie  of  countercurrent  flow. 


Let  us  introduce  tiie  dimensionless  parameters 


1  1  > 

,  ■  I  ,  ■ 


I  I  ■ !  ^'i) 


L,('i  (i^i\ 


-  /.V 

d  '■  1 


Equatiotis  (1)  in  this  case  takcf  the  following  form: 
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E 1  im I na L 1 ng  x'  and  y',  we  obtain  a  system  of  two  differential 
equations  of  the  first  order 
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The  solution  to  this  equation  for  the  case  of  the  cocurrent  flow 
has  the  form 
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where 


<i)  =  I  f  2(1  H  /)P)n  —  (1  — 


q,  =  (1-1-  «)[l  |',(i  -I-  /,|',,ol.  '''* 


III  hhe^e  ovpt'essions  quantities  A'*  and  B*  ai’e  integration  con¬ 
stants,  which  are  determined  from  boundary  conditions. 


In  the  design  of  thermogenerators,  fr'om  our  point  of  view,  as 
the  boundary  conditions  it  is  must  rational  to  accept  the  determining 
values  of  the  maximal  and  mininimal  temperatures  in  the  apparatus. 

Thus  when  m:=  0  we  will  have  x^  =x.T*''  and  y  .  Hence, 

fISKC  r  ’ 

1  1 
■d  *  c  ”(■  i  P  ~  =  P  —  ^^1/-  (7) 

The  heat  fed  to  the  gouierator  and  the  heat  Qx  removed  from 
the  getjeratur,  equal,  respect  i  ve  1  y  ,  to 

<?,  -  -  r, <?.  -  -  7\'n- 

are  determiried  by  means  of  expressions 
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We  fil'd  the  useful  electrical  power  of  the  generator  P  =  0^-Qx, 
eff.  .ency  .jf  the  ‘'.enerator  t]  ==  (G, —  and  average  value  ol' 

1  rie  -utpf.  voltag'  h,  relerped  to  a  sir.gle  thermoe  lement  ,  in  the  form 
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II  solving  the  sysLom  of  equaLions 
’rent  flow,  we  get 


(4)  for  the  case  of  the  cotn- 


J--I-  1 
4  II 


ch 


A\|^2m  (<fi  +  (Ih)  ^4': 


sh  — 


2tl;2  w 

X(cpi  -  q--)  m  ■ 


+ 


(i-«)  H.2"n  r 

1 - ;sl| - J.  OX}!  — 

xl'2  «’>  )  L 

( 1  4~  Aq'2 


y  4'  ^ 


sh - i' 

4‘2  (0 


r  ,  Aq'zrn  (Ti  H-  <p£)  , 

4-  l{  1  ch - 1 - S»1 - 1  >  (’Af 

L  (0  2^2  «■)  }) 


A.((fi  —  q'2) 


OJ 


(11) 


(6a) 


We: 


1  deterniiiie  the  1  nteg*'‘at  ion  constants  from  boundarv  conditions 
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Quantities  A  and  B  liave  the  following  values; 
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Subii  1  tut  i  rig  (1?' 


into  (11),  aid  using  relations  (6)  and  (6a), 
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To  deteriiiine  the  fed  and  v'emoved  quantities  of  heat  Q  and  Q^, 
electrical  power-  P,  efficiency  and  average  value  of  the  output 
voltage  U,  referi’ed  to  a  single  theniioe  lement ,  we  obtain  the  follow¬ 
ing  expressions: 
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It  iR  tjasy  to  Show  that  the  lirruting  ti'aiis  j  t  ion  to  the  case  of 
constant  t  emtc- r'a  t  ure  s  on  each  of  the  tv/o  thecinocoritact  surfaces  (A=0) 
and  the  ahsei;ce  of  thermal  r'esistance  of  the  "therniocontact  pach.ets" 
(^.=  0)  leads  tt:  i'eronulas  unknovni  fr'om  literature  I  1  - 1 J  for  output 
parameters  of  TEG, 

Cuinpaio  son  of  efficiency  of  tiio  cocurt’ent  flow  and  countercur¬ 
rent  flow.  hrt  us  make  a  comparison  with  identical  and 

and  identical  values  of  ji  (i.e,,  at  identical  i.emperatures  of  the 
Ijeai  cai  ricrs  at  the  input  to  the  TEG  and  identical  ourrerit  der'si- 

t  i  e  S  )  . 

L'  L  us  expand  the  cApie.ssiuri  for'  quarititics  6,',  G..",  0,,  f>_,  B- 

iti  series  according  to  powers  of  the  [tararii'-ii-er 
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Usjtig  expression;'  (l6)  for  determining  the  aifferences  in  power- 
arid  efficiency  (with  co jntercurrent  flow  and  cocurrent  flow),  we  get 
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Tilts  the  mode  of  the  countercurrerit  flow  provides  higher  values 
of  ths;  outiiut  power-  and  efficiency  of  the  TEG  than  does  the  mode  of 
tile  coerrt  ci.t  flow. 

Conditiuiis  of  an  ouen  circuit.  When  p  =  0  we  have 

H,-  H/  =-  «'  ((').  «i  =  ^  B:,  --  H  (D). 

Tlie  average  criif  of  the  single  tliei-njoeleri.ent  is  determined  from  ex¬ 
press  i  o  r  i  s 
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1  jf  till  case  uf  a:0,  We  get 
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F.xarri’Tle  of  calculation  of  TEG  paranieLers.  For  the  simplest  case 
G  c  =Gj(C>.  and  k  -ky^  (i.e.  ,  a^O  and  ’u=0),  the  calculation  of  output 
parameters  of  the  TEG  when  1  ,  yf4^f^~0.5,  X  =  0-0.^j  was  conducted; 

f^-0  and  fi=0,l.  Results  of  the  calculation  are  given  on  Figs.  3  and 
4.  It  IS  found  tliat  the  volt-ampere  characteristic  of  the  TEG 
zC  /  a  =  /i(^)  is  close  to  a  straight  line,  and  the  characteristic  of 
the  power  zhP  /  y.I' r  —  is  close  to  a  quadratic  parabola.  Analysis 

of  these  relatioris  shows  that  the  output  parameters  of  the  TEG  (.for 
the  ca.se  ol'  G^c^=Gj.Cj,  arid  k^=kj()  can  be  approximated  by  simple  rela¬ 
tions  having  the  same  form  as  those  known  in  literature  [1-3] 
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c  iT!.jdf  of  maximal  powci-  of  the  TEG  here  corresponds  to  valuc,5 
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It  should  be  rioted  that  the  expressions  obtained  for’  calculation 
of  the  basic  parameters  of  the  TEG  can  be  used  for  solving  other 
problems,  for'  example,  for  an  analysis  of  operation  of  the  TEG  with 
the  earlier'  known  efficiency  of  the  discharge  of  heat  in  the  second 
circuit  arid  so  on . 

Certain  recorrimeridati  ons  on  the  analysis  of  operation  of  the  TEG 
by  taking  into  account  the  dependence  of  thermoelectric  properties 
of  mater'ials  on  temperature.  The  use  of  results  of  the  analysis 
given  above  with  ca i cu lat i on  of  a  specific  design  of  the  TEG  gives 
the  accuracy  accepted  for  engineering  solutions  in  trie  ease  uf  a  rel¬ 
atively  .small  tempei'ature  ciiango  of  tiiermocontact  surfaces  along  the 
ciiaiinel  of  heat  earners.  If’  these  temperatures  vary  greatly,  theri 
tne  absence  of  tiie  calculation  of  a  temperature  dependence  of  the 
therrnoe  1  ec tr' i c  properties  of  the  materials  can  lead  to  considerable 
emrs  in  tiie  calculation.  We  can  avoid  large  ei'rors  by  dividitjg  the 
tiiermogeneratur'  into  sover'al  sections  with  a  relatively  small  (uf  the 
order'  of  10*^)  temperature  cliarige  of  the  tiiermocoritaet  surfaces  and 
perfurmi  rig  a  calculation  for  eacii  section  separ'ately.  Values  of  tem- 
por'utur'os  of  tiie  liut  and  cuJd  heat  carrier-s  at  tlie  beginning  of  tiie 
iirst  section  (wlien  m-O)  can  be  u.sed  as  the  boundar-y  cunditions.  In 
the  cour-se  of  i  i'G  ca  1  c u  1  a 1 1  cin  ,  tiie  heat  being  fed  and  reruuved  ,  tiie 
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ejf^ctrical  power  on  the  section  being  examined,  and  temperatures  of 
heat  carriers  at  the  end  of  the  section  are  determined.  These  tem¬ 
peratures  serve  as  boundary  conditions  for  calculation  oi'  the  second 
section  of  tiie  TEG  with  values  of  properties  of  the  materials  cor¬ 
rected  for  the  new  temperature  mode.  The  calculation  of  the  third 
section  IS  similarly  conducted,  and  so  on. 

Let  us  note  that  for  calculation  of  the  counterfluw  design  of 
the  TEG  over  the  sections,  it  is  necessary  to  introduce  a  change  into 
the  value  of  the  integration  constants  appearing  in  the  equations 
(11).  In  this  case  these  constants  are  determined  from  boundary  con¬ 
ditions  in  the  same  cross  section  of  the  TEG  (when  m=0),  i.e. ,  accord¬ 
ing  to  maximal  values  of  temperatures  of  the  hot  and  cold  heat  car- 
riiu'S.  Here,  instead  of  (12)  we  have 

1  i 

^  =  +  =  //...  // yM«KC  +  1  H - p.  (12) 

2 

Corr'cspond  1  rig  changes  are  introduced  ir.to  the  calculatiori  form- 
u.lns  (I'O-  For  calculatiori  of  the  cocurr'erit  flow  scheme  of  the  TEG 
cV'';r’  the  sections;  all  tlie  analytical  expressions  obtained  earlier 
rcmairi  in  force. 

Gyinbols:  1  -  temperature  of  beat  carrier,  '^K;  t  -  temperature  of 
jurictiori  of  ther  moeiement ,  ;  n  -  specific  ther^mal  emf  of  thermo¬ 

element,  V/deg,  X  -  specific  heat  of  thermoelement,  W/m-deg;  p  - 
.specific  electrical  resistanc  of  thormocl  elment ,  uhms-m:  z  -  param¬ 
eter  of  figur'O  of  merit  of  ther-rnoel  oinent ,  1/deg;  k  -  specific  thermal 
re.s  i  .starice  of  "  ther-mocoritact  packet"  between  the  heat  carrier'  and 
junction  of  tin  thermoelement,  nr^-dog/W;  G  -  mass  flow  rate  of  the 
heat  carrier,  kg/s;  c  -  specil'ic  heat  ;'f  heat  carTier,  J/kg-deg;  h  - 
height  of  thermoelement,  m;  Gp  arjd  .Sj,  -  areas  of  cross  section  of 
branches  of  thermoelernerit ,  n\‘^;  F  -  thormocontact  area,  n/ ;  F^  -  total 
therrriocontact  area  of  TEG,  m'^;  j  -  electrical  curr'ent  density,  A/m'^; 

L  -  average  output  V'  Itage  referred  to  a  single  thermoelernerit,  V;  g  - 
aver'age  emf  oi'  single  tiiermoe  1  ement ,  V;  P  -  output  electrical  i>ower' 
oi'  TEG,  W ;  -  thermal  power-  consumed  by  the  therrnogenerator ,  W; 

Qy  -  thermal  power-  removed  I'rom  therrnogenerator  ,  W;  h  -  effieiency 
oi'  the  ther'iiogerierator- . 

Subscr-!  pts :  r  -  but;  x  -  euld;  p  -  positive  br-ancii;  r'  -  negative 
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branch. 


Superscripts:  Bx  -  input  of  heat 
carrier.  The  output  parameters  of  the 
flow  are  denoted  by  an  asterisk. 


carrier;  etix  -  output  of  heat 
TEG  for  the  case  of  cocurrent 
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